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We present a search for the production of neutral Higgs bosons <f> decaying into t + t~ final states 
in pp collisions at a center-of-mass energy of 1.96 TeV. The data, corresponding to an integrated 
luminosity of approximately 1 fb _1 , were collected by the DO experiment at the Fermilab Tevatron 
Collider. Limits on the production cross section times branching ratio are set. The results are inter- 
preted in the minimal supersymmetric standard model yielding limits that are the most stringent 
to date at hadron colliders. 
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Higgs bosons are an essential ingredient of electroweak 
symmetry breaking in the standard model (SM) . A search 
for Higgs bosons (denoted as </>) decaying to tau leptons 
is of particular interest in models with more than one 
Higgs doublet, where production rates for pp — > <f> — > 
t + t~ can potentially be large enough for observation 
at the Fermilab Tevatron Collider. This situation is re- 
alized in the minimal supersymmetric standard model 
(MSSM) [J, which contains two complex Higgs doublets, 
leading to two neutral CP-even (h, H), one CP-odd (A), 
and a pair of charged {H ± ) Higgs bosons. At tree level, 
the Higgs sector of the MSSM is fully specified by two 
parameters, generally chosen to be Ma, the mass of the 
CP-odd Higgs boson, and tan/3, the ratio of the vac- 
uum expectation values of the two Higgs doublets. De- 
pendence on other MSSM parameters enters through ra- 
diative corrections. At large tan/3, the coupling of the 
neutral Higgs bosons to down-type quarks and charged 
leptons is strongly enhanced, leading to sizable cross sec- 
tions. The Higgs bosons will decay predominantly into 
third generation fermions. 

Searches for neutral MSSM Higgs bosons have been 
conducted at LEP and at the Tevatron 0, H 0. 
These Tevatron searches used between 260 pb _1 and 
350 pb _1 of collider data. In this Letter a search for 
cf> — ► t + t~ with about 1 fb _1 j(| of data is presented. 
At least one of the tau leptons is required to decay lep- 
tonically, leading to final states containing er/ l; fiTh and 
e/z, where Th represents a hadronically decaying tau lep- 
ton. The data were collected at the Tevatron with the DO 
detector between 2002 and 2006 at a pp center-of-mass 
energy y/s — 1.96 TeV. A description of the DO detector 
can be found in Ref . ■ 

Signal and SM background processes are modeled us- 
ing the pythia 6.329 [H Monte Carlo (MC) generator, 
followed by a GEANT-based 0] simulation of the DO detec- 
tor. The signal events are produced with the width of the 
SM Higgs boson. All background processes, apart from 
multijet production and W boson production, are nor- 
malized using cross sections calculated at next-to-leading 
order (NLO) and next-to-NLO (for Z boson and Drell- 
Yan production) based on the CTEQ6.1 [lOj parton dis- 
tribution functions (PDF). 

The normalization and shape of background contribu- 
tions from multijet production, where jets are misidenti- 
fied as leptons, are estimated from the data by using same 
charge e and Th candidate events (er^ channel) or by se- 
lecting background samples by inverting lepton identifi- 
cation criteria {\iTh and efj, channels). These samples are 
normalized to the data at an early stage of the selection 
in a region of phase space dominated by multijet pro- 
duction. The multijet background estimation in the fiTh 
and eTh channels was checked by using an independent 
method to estimate the background: in the ^th channel 
same charge \iTh events were used and in erh channel the 
multijet background was estimated from measurements 



in data of the probability to mis-reconstruct electrons 
from jets. The differences between the estimates were 
used to set the systematic uncertainty on the multijet 
production. The normalization of the background from 
W boson production is obtained from data in a sample 
dominated by W boson + jet events. 

Electrons are selected using their characteristic en- 
ergy deposits, including the transverse and longitudinal 
shower profile in the electromagnetic (EM) calorimeter. 
To reject photons, a reconstructed track is required to 
point to the energy cluster. Further rejection against 
background is achieved by using a likelihood discrimi- 
nant. Muons are selected using reconstructed tracks in 
the central tracking detector in combination with pat- 
terns of hits in the muon detector. Muons are required 
to be isolated in the calorimeter and the tracker ll| . Re- 
construction efficiencies for both leptons are measured in 
data using Z/j* — > p + fi~,e + e~ events. 

A hadronically decaying tau lepton is characterize d by 
a narrow isolated jet with low track multiplicity [12|. 
Three r-types are distinguished: r-type 1 is a single track 
with energy deposited in the hadronic calorimeter (ir^- 
like); r-type 2 is a single track with energy deposited in 
the hadronic and the electromagnetic calorimeters (p^- 
like); r-type 3 is three tracks with an invariant mass be- 
low 1.7 GeV, with energy deposited in the calorimeter. 

A set of neural networks, NN T , one for each r-type, 
has been trained to separate hadronic tau decays from 
jets using Z/j* — > r + r~ MC as signal and multijet data 
as background. The selections on the neural networks 
retain 66% of the Z/"f* — > r + r~ events, while rejecting 
98% of the multijet background. In addition, a neural 
network has been trained with electron MC events as 
background to separate r-type 2 hadronic tau candidates 
from electrons (NN e ). 

The signal is characterized by two leptons, miss- 
ing transverse energy $t and as an enhance- 
ment above the background in the visible mass 



P T2 + Pt) 2 , calculated using the 



four- vectors of the visible tau decay products P T1 2 and of 
the missing momentum Pr = ($T,$x,$yi 0). The com- 
ponents Tfi x and $ y of $t are computed from calorime- 
ter cells and the momentum of muons, and corrected 
for the energy response of electrons, taus and jets. The 
four-vectors of the hadronic taus are calculated using the 
calorimeter for r-types 2 and 3 and the central tracking 
system for r-type 1. 

In the eTh and /j,Th channels, an isolated lepton (e, /i) 
with transverse momentum above 15 GeV and an isolated 
hadronic tau with transverse momentum above 16.5 GeV 
(22 GeV for r-type 3) are required. The pseudorapid- 
ity 1 77 1 is less than 2 for muons and hadronic taus and 
2.5 for electrons. In addition to the background from 
Z/7* — * r + r~ production, a W(— > tv)+ jet event can 
be misidentified as a high-mass di-tau event if the jet 
is misidentified as a hadronic tau decay. The trans- 
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FIG. 1: The distribution of the visible 
mass M v is for (a) p,Th, (b) cr^ and (c) 
e/j, channels. The Higgs boson signal is 
normalized to a cross section of 3 pb. 
The highest bin includes the overflow. 



verse mass, M^/ M = fepfp^ftxO- — cos At/?)] % , is required 
to be less than 40 GeV for the p/Th and 50 GeV for 
the er/j channel. Here, Aip is the azimuthal angle be- 
tween the lepton and T^t- In addition, a selection is 
made in the A(p(e/ [i,, $t) — A(p(t,$t) plane, such that 
Aip(e/n,$ T ) < 3.5 - Atp{T,$ T ) if Atp{T,$ T ) < 2.9 or 
A(p(e/(i,$r) < 0.6 otherwise. This selection removes 
events where the missing transverse energy is in the 
hemisphere opposite to the muon and the tau candidate. 
Due to the larger multijet background in the er^ chan- 
nel the azimuthal angle between the electron and tau, 
Aip(e, r), is required to be greater than 1.6. 

The eT/i channel has a significant background from 
Z/7* — * e + e~ production, where an electron is mis- 
reconstructed as a tau candidate. To remove these 
events, the tau candidates in the er^ channel are re- 
quired to be outside of the region 1.05 < \r)\ < 1.55, 
where there is limited EM calorimeter coverage and are 
required to have less than 90% of their energy deposited 
in the EM calorimeter. Finally, r-type 2 candidates 
are required to have NN e > 0.8, which rejects 92% of 
the Z/-f* — * e + e~ events, while retaining 83% of the 
Z/j* — > r + r~ events. 

We select one muon with pt > 10 GeV and one elec- 
tron with pt > 12 GeV in the ep, channel. Multijet and 
W boson production are suppressed by requiring the 
invariant mass of the electron-muon pair to be above 
20 GeV and $ T + p T + p T > 65 GeV. Background from 
W+jet events can be reduced using the transverse mass 
by requiring that either M T < 10 GeV or M% < 10 GeV. 
Furthermore, the minimum angle between the leptons 
and the $t vector, min{Atp(e, $t) , Aip(p, $t)], has to be 
smaller than 0.3. Contributions from tt background are 
suppressed by rejecting events where the scalar sum of 
the transverse momenta of all jets in the event is greater 
than 70 GeV. 

The number of events observed in the data and ex- 
pected from the various SM processes show good agree- 
ment (Table . The number of background and signal 
events depend on numerous measurements that introduce 
a systematic uncertainty: integrated luminosity (6.1%), 
trigger efficiency (3%-4%), lepton identification and re- 
construction efficiencies (2%-10%), jet and tau energy 
calibration (2%-3%), PDF uncertainty (4%), the uncer- 



TABLE I: Expected number of events for backgrounds, num- 
ber of events observed in the data and efficiency for a signal 
with M<f, — 160 GeV for the three channels. The uncertainties 
are statistical. 



Channel 


er h 




e/i 


Z/-f* — > T + T~ 


581 ±5 


1130 ±7 


212±3 


Multijet 


332 ± 20 


86 ±4 


29 ±1 


W — * ev, fiv, tv 


42 ±5 


32 ±4 


9±2 


Z/j* — > e + e~,/i + /i~ 


31 ±2 


19 ±1 


12 ±1 


Diboson + it 


3.0±0.1 


7.0 ±0.4 


6.1 ±0.1 


Total expected 


989 ± 23 


1274 ± 9 


269 ±3 


Data 


1034 


1231 


274 


Efficiency (%) 


1.04 ±0.03 


1.46 ±0.04 


0.57 ±0.03 



tainty on the Z/j* production cross section (5%), nor- 
malization of the W boson background (6%-15%), and 
modeling of multijet background (4%-40%). All except 
the last one are correlated among the three final states. 
Most of the uncertainties affect only the overall accep- 
tance for the signal and backgrounds. However, uncer- 
tainties on the energy scale and electron trigger efficien- 
cies modify the shape of the visible mass distribution 
(Fig. [T]) . These uncertainties are therefore parameterized 
as a function of M w - ls . 

We extract upper limits on the production cross sec- 
tion times branching ratio as a function of Higgs boson 
mass M,),. In order to maximize the sensitivity (median 
expected limit), the event samples of the er^ and p/Th 
channels are separated by r-type to exploit the differ- 
ent signal-to-background ratios. Furthermore the differ- 
ences in shape between signal and background are ex- 
ploited by using the full M vis spectrum in the limit cal- 
culation (Fig. []}. The limits are calculated by utiliz- 
ing a likelihood-fitter [l3[ that uses a log-likelihood ratio 
test statistic method. The confidence level, CL S , is de- 
fined as CL S = CL s+ b/CLb, where CL s+ b and CLb are 
the confidence levels in the signal-plus-background and 
background-only hypotheses respectively. The expected 
and observed limits are calculated by scaling the signal 
until 1 — CL S reaches 0.95. The resulting cross section 
limits are shown in Fig. [21 The difference between the 
observed and expected limits at high masses is slightly 
above two standard deviations. It is mainly caused by a 
data excess in the \itk channel above M v - m of 160 GeV. A 
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large number of kinematic distributions were studied for 
this sample and the data are consistent with both back- 
ground and signal shapes. Due to the M vis resolution 
these events affect the limit over a wide range of masses. 
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FIG. 3: Ratio of expected cross section limits using a Higgs 
boson with non-SM width to those calculated with a Higgs 
boson with SM width, as a function of F^/M^. 



FIG. 2: Expected and observed 95% CL upper limits on the 
cross section times branching ratio for </> — > t + t~ production 
as a function of M$ assuming the SM width of the Higgs 
boson. The ±1,2 standard deviation bands on the expected 
limit are also shown. 



currently disfavored [161 ]. The production cross sections, 
widths, and branching ratios for the Higgs bosons are 
calculated over the mass range from 90 GeV to 300 GeV 
using the feynhiggs program 14[. In these scenarios 



The limits in Fig. [2] assume a Higgs boson with SM 
width, which is negligible compared to the experimental 
resolution on M v - m . In models such as the MSSM the 
Higgs boson width can become substantially larger than 
the value in the SM. This was simulated by multiplying 
a relativistic Breit-Wigner (BW) function with the cross 
section from FEYNHIGGS [lj] for masses M > 80 GeV 
to obtain the differential cross section for a wide Higgs 
boson as a function of mass: 

= <r(M, tan A I> = 0) x BW{M,M^). (1) 

This differential cross section was used to build a signal 
template of the M vis distribution for a Higgs boson of 
mass M$ and width L^. The limit calculation procedure 
was then repeated with templates corresponding to vari- 
ous values of T^. The ratio of the expected cross section 
limit for a wide Higgs boson to the limit for a Higgs bo- 
son with SM width as a function of L^ /M$ is shown in 
Fig. [31 This result can be used to correct the cross sec- 
tion limit for a Higgs boson with SM width (Fig. [2]) for 
a non SM width in a model independent way. 

In the MSSM, the masses and couplings of the Higgs 
bosons depend, in addition to tan/3 and Ma, on the 
MSSM parameters through radiative corrections. In a 
constrained model, where unification of the SU(2) and 
U(l) gaugino masses is assumed, the most relevant pa- 
rameters are the mixing parameter X t , the Higgs mass 
parameter /i, the gaugino mass term M%, the gluino 
mass m g , and a common scalar mass MsuSY- Limits 
on tan/3 as a function of Ma are derived for two sce- 
narios assumin g a CP-conserving Higgs sector 113: the 
TO max scenar i hj| anc j th e no-mixing scenario [19( with 
fi = +0.2 TcV. The /i < case is not considered as it is 



m™ ax , n = +200 GeV No-mixing, \i = +200 GeV 




M A (GeV) M A (GeV) 



FIG. 4: Region in the (Ma, tan 0) plane that is excluded 
at 95% CL for the m™ ax and the no-mixing scenario 
(m t = 172.6 GeV Q3). Also shown is the excluded region 
from LEP 0]. 

Ta/Ma < 0.1 for M A < 200 GeV. The effect of the 
Higgs boson width is therefore small. For large tan (J, 
the A boson is nearly degenerate in mass with either the 
h or the H boson, and their production cross sections 
(gg —t<j),bb—t <j>) are added. 

Fig. [4] shows the results interpreted in the MSSM sce- 
narios considered in the Letter. We reach a sensitivity of 
around tan/3 = 50 for M A below 180 GeV. The result 
represents the most stringent limit on the production of 
neutral MSSM Higgs bosons at hadron colliders. 
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